A PSEUDO RANDOM NUMBER GENERATOR FROM CHAOS 

NABARUN MONDAL AND PARTHA P. GHOSH 

Abstract. A random number generator is proposed based on a theorem about 
existence of chaos in fixed point iteration of a;„ — cot^{xn-i)- Digital computer 
simulation of this function iteration exhibits random behavior. A method is 
proposed to extract random bytes from this simulation. Diehard and NIST test 
suite for randomness detection is run on this bytes, and it is found to pass all 
the tests in the suite. Thus, this method qualifies even for cryptographic quality 
C^ random number generation. 
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1. Introduction 



— 1 Fast generation of high quahty random numbers using simple arithmetic meth- 

CsJ ods are elusive. The most elusive aspect about generating random numbers is that 

^_^ one can never be sure, if the produced numbers are really random. That is exactly 

^ why a background theory is needed for it, and a set of standardized statistical tests 

^ must certify the numbers as random. 

. Robert R. Coveyou suggested : "The generation of random numbers is too im- 

Q portant to be left to chance. ". Almost same theme is reiterated in the immortal 

'""' words of Knuth [4J "one should not expect arbitrary algorithms to produce random 

t^^ numbers, a theory should be involved" . The great Von Neumann asserted ran- 

^ domness lies beyond the realm of arithmetic: "Any one who considers arithmetical 

^r^ methods of producing random digits is, of course, in a state of sin. " As it turned 

I^ out to be. Von Neumann might not had been exactly correct, and that is what the 

. authors want to present in the current paper. 
^ Until recently it was assumed that randomness comes out of complexity. How- 

psj ever, it can also be an effect of chaotic dynamics [fj- Properties of a chaotic 

^ system would be discussed in the section [2J Many people misleadingly talk about 

J> as if there are various sources of randomness. According to some authors [7] , 

KA the fundamental source of randomness, if at all any, is Heisenberg Uncertainty 

H Principle. 
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A practical viewpoint is : Randomness occurs to the extent that something can 
not be predicted. Randomness is a matter of degree, and lies in the eyes of the 
beholder. Poincare pointed out that the classic random outcome of a die throwing 
or a flipping coin, comes from sensitive dependence to the initial condition. A 
small perturbation causes a large difference in the final outcome, thereby making 
prediction difficult. This sensitive dependence is a hallmark of Chaos [lOj . 

Chaotic trajectories even look random, and they pass many classic "tests" of 
randomness. This in fact generates the principle of equivalence between chaotic 
and random systems, as discussed in [1]. 

Furthermore, chaotic systems might arise in very simple forms of iterative maps 



O [in] (definition 2.1). Two such simple looking systems are discussed in the 
section [2| both exhibiting complex behavior. They form the operating principle 
(or the theory) behind the proposed random number generator is presented. 

However, chaotic systems trajectories sometimes can approach the subset of the 
state space called attractors, thereby drastically reducing the perceived random- 
ness. This is a natural outcome of running the simulation in digital computer, 
where the precision of the state point is limited to fixed number of bits of infor- 
mation, without arbitrary precision, the simulated chaotic system would get into 
a cycle. 

The practical principle of garnering randomness out of chaotic systems would 
be then, to ensure, that the system trajectory does not end up being in an attractor. 
Practical generation of random digits (bytes) are the topic of discussion in section 

m 

Finally, the result of two of the industry standard statistical test suites are 
presented in section |4j It was found that the generated random byte stream passes 
all these industry standard tests, thereby making the generated numbers suitable 
even for cryptography. 



2. Theory of Operation 

In this section we would discuss the operating principle of the proposed random 
number generator, and prove that, the iterative map is Devany Chaotic (definition 



2.2 [3] ). Some of the prerequisite terms and definitions can be found in Appendix 

w 

We start with the idea of a fixed point iteration, or a single dimensional discrete 
dynamical system. 



Let there be a function / : X — )■ X where X is a metric space (definition A.2). 
Lets start with xq € X, and let Xi = /(xq). Let X2 = f{xi). 

Definition 2.1. Iterative Function. 

The general form :- 
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is called iterative function system, or a fixed point iteration, which is a type of 
single dimensional discrete dynamical system, or map. 

Many iterated systems exhibit what is known as "chaotic behavior". "Chaos" 
however remains a tricky thing to define. In fact, it is much easier to hst properties 
that a system described as "chaotic" has rather than to give a precise definition of 
chaos |3|. 

Definition 2.2. Chaotic Systems. 



A chaotic fixed point iteration (dynamical system 2.1 ) is generally characterised 
by:- 



(1) Having a dense (definition A. 6) collection of points with (definition A. 3) 
periodic orbits. 

(2) Being sensitive to the initial condition of the system (so that initially nearby 
points can evolve quickly into very different states), a property sometimes 
known as the butterfly effect. 



(3) Being topologically transitive (definition A.l). 



2.1. Discontinuous Maps. The map we would be discussing here is a; ^ cot^{x). 
This is a discontinuous map. Unfortunately, while literature is filled with material 
on continuous maps, very few are really available on discontinuous maps [2]. Re- 
ciprocal of any continuous map becomes discontinuous on the set of zeros of the 
original map. 

For example the map x <— l/x^ is a reciprocal map of x^, while the zero of the 
original a; = 0, became a point of discontinuity. 

In general there is nothing interesting about these type of reciprocal maps, as 
this specific one has a fixed point x* = ±1, and while iterated would immediately 
converge to one of them. 

But, these maps behave very differently when put into a circle domain, that is 
in modular form. 

Such a modular form for the map x ^ 1/x^ is:- 

Xn = I -^ — ) rnod 1 ; x G (0, 1) 

For some maps which are inherently periodic, this contraption becomes a natural 
choice to investigate their behaviour. 

For example we can define the reciprocal map x -^ l/cos{Xn~i), whose domain 
is :- 

D = R\{{2n + l)7v/2} ; n e Z 

But what happens in real is modular arithmetic with period tt. Hence, any output 
x„ > vr or Xn < should be treated modulo vr, due to periodicity of cos{x). In this 
case, the modified system have the domain changed from D to :- 

D^ = [0,7i/2)U{n/2,n] 
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The output range should also be changed accordingly to :- 

Xn = 7 r rnod n 

COS[Xn~l) 

This works in a domain where 0, tt are actually the same point. Notice the un- 
necessary use of the period vr in the definition. Therefore, we can generalise this 
for any periodic function with a period T and we can go a step beyond to always 
normalise the domain and range as the same interval: a.k.a the unit circle :- 

/ = (0,1) 

by introducing the following definition of normal niaps:- 

Definition 2.3. Normalised Maps. 

Let g{x) be any function with a period T . Then:- 

iV»W = ^^^^^;x.(o.i) (2.1) 

is called the normalised rational map form of g{x). 

With this definition is mind, lets see how a modulus operation applies on a 
linear map y{x) = mx + C with mod T. We note that obviously when yn = T , 
yn would become zero, and would generate discontinuity, but while yn < T, the 
derivative remains y' = m. The modulus operation does not change the derivative 
of the underlying function, just makes it undefined at periodic points nT ; n & Tj. 



Comparing equation (2.1) of definition 2.3 with the sawtooth map equation (2.2) 



suggests that the map y{x) would become a sawtooth type map, which is known 
to be of chaotic type, given the slope is greater than 1, that is |r7i| > 1. 

/(x) = m(l - x) mod 1 ; m G M (2.2) 

In case of a general function g{x), only the slope would vary, which would become 
a function of x itself as in :- 



// 



m[x) = g [X) 



A normalized function (definition 2.3) diverging to infinity (having singulari- 



ties) would still have the derivative of its non modular analog, but would have 
introduced countable discontinuities due to modular operation. 

Image ([I| demonstrate the effect modularisation has on the function 1/x^. Note 
the discontinuities appearing on the side of x — )■ 0, they are all chopped off to 
(0, 1), and as we are mapping a potentially infinite length with unit interval, we 
would need countable amount of them. 

The fixed point x* = 1 can not be reached from anywhere in the domain now. 
If an iteration value Xn gets close to the value 1, in the next iteration it would be 
pushed back to the left side of the map, with x„+i ~ 0. Which means it would go 
into the zone of the chaos, where it will be eventually pushed back to any point 
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other than 1. This dynamics is what essentially generates the chaotic behaviour in 
the system. These type of systems have tremendous sensitivity to initial condition. 




Figure 1. Function : l/x^ mod 1, showing the sawtooths 

In fact, the modularisation makes many trivial map act in a non trivial way, for 
example the always 3-cycle map of a; ^ 1/(1 — x) under the normal form 

mod 1 ; X G (0, 1) 



Xr. 



1 — X 



becomes just as chaotic as depicted in Image [Tj 

In this regard, it is curious to note how modulus operation can generate chaotic 
behaviour. Modulus operation is a very special case of the general discontinuous 
maps C studied in [2], which we define next. 

Definition 2.4. Sharkovosky Chua Type Map [2j •' f ^ C. 

Let f : I ^f I with [a, h] he map such that:- 

(1) f is continuous everywhere except points :- 

D = {zi,Z2,...,Zr} ; r e Z 

(2) f is monotonic in the interval {zi,Zi+i) for i = {0, l,...,r}, with zq = 
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(3) The limits 



and 



lim /(x) ^ lim /(x) 

x^Zi~0 x—^Zi+0 



lim = f(x) G {a, b}\/i 

x—^Zi±0 



(4) It is expansive |/'(x)| = / > 1. For every point x G / \ {zi} there exists an 
interval U^ containing the point x such that:- 

d{f{U)) > ld{U) 

for every interval U G Ux where d{V) is the length of the interval V. 
then, f E C is called a Sharkovosky Chua Type Map. 

Note that the original papers P] definition does admit countable number of 
infinities, but that admission does not change any property of the map with respect 
to chaotic behaviour, in fact augments it. 

Now, we would introduce a special class of maps with more general proper- 
ties then Sharkovosky Chua type maps, these maps allow a specific type of non 
expansive intervals in the function domain. 

Definition 2.5. Generalized Sharkovosky Chua Type Map : f E Q. 

A function f is said to be of Generalised Sharkovosky Ghua type f E Q, iff the 
function has the properties from (2,3) and modified (1, 4) with the following:- 

(1) Gountable number of discontinuity. 

(2) Every element Xg E Ig such that \f'{xs)\ < 1, has an element in forward 
orbit Xg such that \f'{xg)\ > 1. 

(3) It is expansive |/'(x)| = / > 1. For every point x E I \ {{zi} U {Is}) there 
exists an interval U^ containing the point x such that:- 

d{f{U)) > ld{U) 

for every interval U C Ux where d{y) is the length of the interval V . 

Lemma 2.1. Reduction to Sharkovosky Ghua Type Map. 



A Generalised Sharkovosky Ghua Type map (definition 2.5 ) can be contracted 
into a Sharkovosky Chua Type map (definition 2.4^- 



Proof. The reduction is straight forward. 

Informally, we contract the orbits from an expansive region to a non expansive 
region, to another expansive region. 

As the intervals {Ig} forwards all the points to the points Xg such that \f{xg)\ > 
1, we can contract those iterative steps (the orbit from Xg> — ?■ x^ — > Xg ) into 
a single step, and eliminate the whole intervals {Ig} where |/(xs)| < 1, thereby 
contracting the intervals, and no point Xg remains where \f'{xs)\ < 1. 

D 
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Theorem 2.1. Sharkovosky Chua Type Maps are Chaotic. 

Sharkovosky Chua Type Maps are chaotic, so are the generalised Sharkovosky 
Chua Type maps. 

Proof. This has been discussed in j2]. As the generahsed Sharkovosky Chua maps 
can be contracted into a Sharkovosky Chua map, they are chaotic too. D 

Now we give one example of a normahsed map that is not in Q. Curiously , the 
map l/sin^{x) ^ Q, because the non expansive zone does not push the orbit back 



to the chaotic zone, as required by definition (2.5) property (3). As the fixed point 



(about 0.37) lies in the non expansive zone, the map converges, and is not chaotic 
at all. 

This can nicely be explained with the Image |2j see at the zone where the function 
has zero, does not have the expansive property, the zone is stable. 
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Figure 2. Normal form of [ l/sm^(x) - x ] showing the fixed point zone 

What maps are of Generalised Sharkovosky Chua type / G ^ ? Take a polyno- 
mial map with singularity, normalise it, it probably goes in Q. Examples are l/x^ 
map, 1/(1 — x) map. The trigonometric reciprocal maps l/cos{x) G Q, so is the 
map l/cos^(x) The map tan{x) G Q and cot{x) G Q. 

Finally, the map cot'^{x) G Q, which is discussed in the next subsection. 
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2.2. Reciprocal Cot Squared Map. In this subsection we discuss tlie dynamics 
of the cot~'^{x) map. This is formally represented by equation (2.3) 



Xr, 



COs{Xn-l) 



(2.3) 



Cof^i^Xn-l) \Sin[Xn-l] 

This map can is presented in the semi normal form, that is the domain as (0, 1.0) 
but the range is in M in the figure ([3]). In the figure dsl), we did not present the 
singular points {0, 1} and draw only the range (0.1, 0.9). 
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Figure 3. The function cot^(x7r) 

It is apparent that due to symmetry, to normalise the function it can be taken 
as in 7r/2 instead of vr, as if we normalise with period tt, the normalised function 
N^ would be Nc{x) = N^{l-x). 

The function normalisation can be seen in figure dll). 

Theorem 2.2. Iteration of Xn+i = cot'^{xn) is chaotic. 

The iteration of Xn+i = cot^{Xn) by being a generalised Sharkovosky Chua type 
map \2.1\ is chaotic. 



Proof. The normalised function (definition 2.3) form for cot [x) with period 7r/2 



can be written by equation (2.4) 
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Figure 4. Normal form of cot'^{x) from equation (2.4) 
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Figure 5. The zone of chaos in Nc (equation 2.4) 
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Figure 6. The Normalised Derivative cot'^{x) (2.5) 



cot^ — I 



TT 



N,{x) = ^l^cot'i^^j mod-) (2.4) 

It is apparent from the figure ^ that for values of x„ — )■ 1 the next value a;„+i 
would be in the chaotic zone, depicted in figure ([s]). 
The derivative is :- 

d[cot'^{x)] _ 2cos{x) 



dx 



sin^(x) 



(2.5) 



We note that the derivative |/'(x)| < 1 only when x„ > 0.7 however, at that 
point the next Xn+i would push it between (0.1,0.3) which is a massively chaotic 
zone as depicted by figure (IS]). Hence, every non expansive point maps back into 
the zone of chaos, which is also the expansive zone with |/'(a:)| > 1 . 

Hence equation (2.4) satisfies condition of generalised Sharkovosky Chua type 
map, and hence is chaotic. D 

We now present the theorem of principle of equivalence between chaotic and 
random behavior, which would be used to apply the random behaviour of the 
iterate cot'^{x). 

In a very insightful paper pj it has been proven that:- chaotic and random 
systems are observationally indistinguishable. If that is the case, then, one can 
replace a random system by an equivalent chaotic system, and vice versa, as has 
been argued in \T\. 

Theorem 2.3. Chaotic systems and Random Systems are Observation- 
ally Equivalent [1\. 
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A chaotic system can be observationally replaced with a random system, and vice 
versa. 



3. Experimental Validations 



As the theorem 2.3 suggests, by virtue of being chaotic, the iteration of theorem 
(2.2) could serve as a randomness generator. However, there are cautionary words 
involved because arbitrary precision is needed while doing so. 

When the digits are of finite precision, say 2 digits before and after decimal, 
the system can not exhibit an arbitrary long period, as the total number of states 
possible are lO^"*"^ = 10^. In general if the number of digits in base b are of N , 
then the maximum possible number of states 15*1 becomes:- 15*1 = b^ which is still 
finite. However, as we keep on increasing A^, IS"! increases exponentially, and the 
chance of repeat decreases. 

The perl code for computing the iteration of x = cot'^{x) is presented here. 



1 


sub cot 




2 


■c 


^^m 


3 


my $x = $_ [0] ; 


^^M 


4 


cos($x) /sin ($x) 


■ 


5 


} 




6 


while ( 1 ) 


^^1 


7 


i 


^^1 


8 


printf "7..18f\n", 


^H 


9 




$x = cot($x) ** 2 
} 


' 



The above code generates x G IR+, which is not well suited for visualizing the 
generated randomness. A chaotic dense orbit can only be visualized in case of a 
bounded domain. To move the resulting value in the bounded (0, 1) domain the 
following code can be used:- 



while 


( 1 ) 




{ 


$x = cot($x) ** 2 ; 


J 




my $fract = $x - int ( $x 


^ fl 




printf "7..18f\n", $fract 


■ 


> 







This code generates x G (0,1). The result of the iteration is dense (definition 
A. 6) as depicted in the figure [tI and a close-up view of the iteration is shown in 



the figure [8j 

However, all industry standard statistical tests run on binary byte streams of 
data. Therefore to test for randomness one has to transform this sort of random 
"real" data points into a stream of random bits. 
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Figure 7. Iteration of x„+i = cot^{xn) (fractional part) 




Figure 8. Iteration of x„,+i 
(fractional part) 



COt^iXr 



Showing a small Range 



Hence, we needed to convert this apparently random iteration values into stream 
of bytes so that randomness test suites like diehard/NIST can be used to verify 
the randomness. This is discussed in the next subsection. 

3.1. Digital Representation of Real. The IEEE floating point specification 
mandates that the double representation would be of 64 bits or 8 bytes. If the 
representation starts with bit number as 6o and ends with bit number 63 as fegs ■, 
then, the &62--^52 is exponent of the number, while fegs is the sign bit. In byte wise 
notation it is BjBq...BiBq . 

So, if the function x„+i = cot'^{xn) is really chaos generating, then, if we can 
gather the mantissa : 651... 6o of the x„ , the resulting bitstream should behave as 
random. However, it wold not be optimal for extracting 651. ..60 , as they are not 
in the byte boundaries. A better, faster option would be discarding 2 of the most 
significant bytes, that is BjBq and taking the other bytes B5B4...B0 as random 
bytes. 
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The Python language uses hex{) call to show the floating point number :- 

[sign] [Ox] integer [. fraction] [p exponent] 

Hence, the following python code demonstrates the exact random bytes as in the 
hexadecimal string. 



1 


im 


'P 


ort 


math 








2 


im 


'P 


ort 


time 








3 


X 


= 


time.timeO #Seed 


it , 


go to 


Chaos 


4 


X 


= 


X 


- math . floor ( x 


) 






5 


wh 


i 


le 


True : 








6 






y 


= math. tan (x) 








7 
8 
9 
.0 
.1 
2 






y 

X 

s 
a 
b 
s 


= y * y 

= l.O/y 

= X . hex () 

= s. split (".") 

= a[l] . split ("p" 

= b[0] 


') 






-3 


1 




print s 









The resulting output is shown below:- 



1 49c07b03923a2 

2 2201d34423ec8 

3 47527dadaecb0 



private double number ; 
public DoublelmplO 
{ 

number = System . nanoTime () ; 
} 

public double getNumberO 
{ 

double r = Math . tan ( number ) ; 

number = l/(r*r) ; 

return number ; 
} 

public byte [] getRandomBytes () 
{ 

long 1 = Double . doubleToRawLongBits (getNumber ()) ; 

byte [] bytes = toByta_Str ipSignByte (1 ) ; 

return bytes ; 
} 

/* Return bytes from B_5 to B_0 here */ 
public static byte [] toByta_Str ipSignByte ( long data) ■[ 

return new byte [] { 

(byte)((data >> 40) & Oxff), //B_5 
(byte) ((data >> 32) & Oxff), //B_4 
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(byte) ((data >> 24) & Oxff), //B_3 

(byte) ((data >> 16) & Oxff), //B_2 

(byte)((data >> 8) & Oxff), //B_l 

(byte)((data >> 0) & Oxff) //B_0 




The above code is platform independent, and authors also wrote a faster C imple- 
mentation. 

4. Results of The Statistical Testing 

For statistical testing suites a binary file is to be passed to the test suite, which 
is then read as a stream of random bits. 

We called the getRandomBytes function (which generates 6 random byte each 
time it gets called) to generate medium (min:600 Mega Bytes) to long ( Max: 6 
Giga Bytes) binary files. 

On a 2.66 GHZ Core 2 Duo Macbook pro, with OS X Lion, using the C langauge 
variant , the 60 Megabyte file takes 2 sec, for 600 MB 19 sec, and 6 Giga Bytes of 
random bytes generation takes 290 sees. 

This file, was passed to the diehard family [5j of test suite. The findings are 
summarised in the table [I]. 

However, diehard is an old generation of test suite, and NIST suite [6] is the 
current industry norm for randomness testing. This generator passes all the NIST 
suite |;6j tests with ease. The file size of 6 gigabyte was passed to the NIST suite, 
with minimum bitstream length of 10000 bits, with minimum number of bitstream 
1000, to a maximum of 10000. The resulting p-values are aggregated over all the 
bitstreams. 

The results are summarized in the table [2]. 



5. Usage Advantages 



Due to 



The premise of the operation of the cot^ generator is the theorem 2.2 
the nature of the iteration (chaotic), the following properties are true :- 

(1) Hard to Predict. 

The seed of the generator is xq, which can be an arbitrary precision 
floating point value. For the orbits of the iteration with Xq = 1.0000000 
and with Xq = 1.0000001 the resulting byte stream start differing from 
the 3rd byte. The choice of the default initial value xq as current systems 
time in nano second, makes it impossible to track and predict the value 
for any iteration, unless absolutely precise time synchronisation is achieved 
between two systems. 

Due to machines architectural and implementation differences, the same 
code produces two completely different byte streams in C++ and Java. 
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Test Name 


Resulting p-value[s 


End Result 


Birthday Spacing 


0.948788 


PASSED 


Overlapping 5-Permutation 


0.874604,0.993319 


PASSED 


Binary Rank Test (31X31) Matrix 


0.724738 


PASSED 


Binary Rank Test (32X32) Matrix 


0.853068 


PASSED 


Binary Rank Test (6X8) Matrix 


0.028129 


PASSED 


The Bit Stream Test 


min 0.0545, max 0.8906 


PASSED 


OPSO 


min 0.0373, max 0.9579 


PASSED 


OQSO 


min 0.0500, max 0.9688 


PASSED 


DNA 


min 0.0359, max 0.9986 


PASSED 


Count The I's On Byte Stream 


0.223283, 0.165772 


PASSED 


Count The I's On Specific Bytes 


min 0.022031, max 0.967922 


PASSED 


Parking Lot 


0.343578 


PASSED 


Minimum Distance Test 


0.097679 


PASSED 


3D Sphere Test 


0.885174 


PASSED 


Squeeze Test 


0.805752 


PASSED 


Overlapping Sums Test 


0.531035 


PASSED 


Runs Test Up 


0.297245, 0.766326 


PASSED 


Runs Test Down 


0.315770, 0.857272 


PASSED 


Craps Test 


wins:0.360479,throws/game:0.207756 


PASSED 



Table 2. Results of the NIST Suite 



Test Name 


Resulting p- value s 


End Result 


Frequency 


0.379045 


PASSED 


Block Frequency 


0.963497 


PASSED 


Cumulative Sums (2 tests) 


0.660844,0.895204 


PASSED 


Runs 


0.013760 


PASSED 


Longest Run 


0.928150 


PASSED 


Rank 


0.007694 


PASSED 


FFT 


0.081013 


PASSED 


NonOverlappingTemplate (multiple) 


- 


PASSED 


Universal 


0.678686 


PASSED 


Approximate Entropy 


0.048404 


PASSED 


RandomExcursions (multiple) 


- 


PASSED 


RandomExcursions Variant (multiple) 


- 


PASSED 


Serial (2 tests) 


0.019966, 0.175884 


PASSED 


LinearComplexity 


0.679514 


PASSED 
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Unless one knows the exact architecture of the machine and environment 
(language of implementation) on which the algorithm is running, it would 
be hard to predict the outcomes. 

(2) Conjectured to be Unique. 

For sufficient precision , the value of the Xn will almost never repeat. 
The machine simulation shows this to be true for default precision, up to 
millions of iterations. Given the whole system can be made using arbitrary 
precision floating point arithmetic libraries of choice, then the values would 
be almost surely unique. 

(3) Fast Algorithm. 

The Java implementation for this random generator is faster than 

Java, security . SecureRandom 

implementation. To generate a 600 mega byte of random bytes, using 6 
bytes of buffer, current Java implementation takes 150 seconds on a SUN 
Ultra-SPARC machine, running Java 7 on SOLARIS, while SecureRandom 
takes more than 230 Seconds. 

6. Summary and Future Works 

In this paper we have demonstrated a computationally easy and fast method of 
generating very good random numbers for all practical purposes, including cryp- 
tography. Arbitrary large random chunk of bytes can be generated if we take 
an arbitrary precision Real Number implementation like BigDecimal in Java, al- 
beit the number generation would get slower. On the flip side, as the principle is 
well established, we can guarantee arbitrarily large period for the random number 
generator. 

Appendix A. Definitions of The Theory Section. 

Definition A.l. Fixed Point of a function. 

For a function f : X ^ X , x* is said to be a fixed point, iff' f{x*) = x* . 

Definition A. 2. Metric Space. 

A metric space is an ordered pair (M, d) where M is a set and d is a metric on 
M , i.e., a function:- 

d: M X M ^R 

such that for x,y,z & M , the following holds:- 

(1) dix,y)>0 

(2) d{x,y) = iffx = y . 

(3) d{x, y) = d{y, x) 

(4) d{x, z) < d{x, y) + d{y, z) 

The function 'd' is also called "distance function" or simply "distance". 
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Definition A. 3. Orbit. 

Let / : X — )■ X be a function. The sequence O = {xq,Xi,X2,X3, ...} where 

Xn+i = f{xn) ; Xn e X ; n > 
is called an orbit (more precisely 'forward orbit') of f . 

Definition A. 4. Period of An Orbit. 

The length of the orbit (definition A.3 ) is called the period of the orbit. Hence, 
for an orbit O the period 'p ' is:- 

p= \o\ 

f is said to have a 'closed' or 'periodic' orbit O if \0\ ^ oo or the period is not 
infinity. 

Theorem A.l. Convergence Criterion for a Fixed Point Iteration (Ba- 
nach Fixed Point). 

Iteration of a function / : X — )■ X Xn+i = f{xn) would converge to the fixed 
point X* iff:- 

|/(x)-/(x*)|<|x-x*| 
where x, a;* G A C X where A is called the basin of attraction. 

That means, in short if x* be the fixed point, and |/'(a;)| > 1 in the neighbour- 
hood of the X* then the fixed point iteration won't converge at the fixed point 

X . 

Definition A. 5. Topological Space. 

Let the empty set be written as : 0. Let 2^ denotes the power set, i.e. the set 
of all subsets of X . A topological space is a set X together with t C 2^ satisfying 
the following axioms:- 

(1) (/} er andX er , 

(2) T is closed under arbitrary union, 

(3) T is closed under finite intersection. 
The set r is called a topology of X. 

Definition A. 6. Dense Set. 

Let A be a subset of a topological space X . A is dense in X for any point x E X, 
if any neighborhood of x contains at least one point from A. 

The real numbers M with the usual topology have the rational numbers Q as a 
countable dense subset. 

Definition A. 7. Topological Transitivity 

A function f : X -^ X is topologically transitive ,if, given any every pair of non 
empty open sets A,B (Z X , there is some positive integer n such that 

f^{A)r}B^^. 

where /" means n'th iterate of f . 
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